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Abstract 

Reformulation of Dirac equation in terms of real quadratic division 
algebra of quaternions is given. Similar equations with different mass 
term are identified as suitable for description of free propagating quark 
motion. The complete orthonormal set of the positive-energy plane 
wave solutions is presented. Therefore, Zitterbewegung phenomenon is 
absent in this formulation. The probability current is proportional to 
the momentum, as in standard Schrodinger wave mechanics. 

PACS: ff.10 Eq 

1 Introduction 

In this paper I will discuss an old, well-known and apparently "insignificant" 
problem that lies in the foundation of the relativistic quantum mechanics or, 
more precisely, in the foundation of the quantum electrodynamics. A story 
starts with the celebrated paper of P. A.M. Dirac that describes relativistic 
motion of a free propagating electron, and had tremendous success of natural 
explanation of spin, correct non-relativistic limit, correct coupling with the 
external magnetic field, correct gyromagnetic ratio and, finally, prediction 
of positron. 
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The mathematical formalism used requires presence of the negative ki- 
netic energy states in order to obtain the complete orthonormal set of the 
linearly independent fundamental solutions of the suggested equation. In- 
deed, such states do not make sense from the physical point of view. 

Since the Dirac equation is written in the Hamiltonian form, it allows us 
to work in the Heisenberg representation and determine directly whether or 
not a given observable is a constant of motion. Then we find, for example, 
that the momentum is a conserved quantity as it should be. However, the 
orbital momentum, as well as the spin, are not conserved separately and 
only the sum of them is a constant of motion. Conventionally, the spin is 
associated with the internal degree of freedom of the electron and therefore 
apparently has nothing to do with isotropy of the space-time continuum 
(indeed, if we assign to the quantum mechanical space-time point internal 
algebraic structure, then this will be rather naturally expected result). Even 
more surprising result [2. is obtained if we consider velocity x in the Dirac 
formulation. Instantaneous group velocity of the electron has only values 
±c in spite of the non-zero rest mass of electron. In addition, velocity of a 
free moving electron is not a constant of motion. 

An analytic solution for the coordinate operator of a free propagating 
electron was found by E. Schrodinger |3J. It turns out that in addition to 
the uniform rectilinear motion consistent with the classical electrodynamics, 
the Dirac electron executes oscillatory motion, which E. Schrodinger called 
Zitterbewegung. Let us recall that entire non-relativistic quantum mechanics 
was raised in order to explain the absence of radiation during the oscillatory 
motion of the electron bounded by the electric potential of the nucleus. 
Therefore, the Dirac theory of electron contains a definite prediction that the 
free moving electron will loose all a kinetic energy through electromagnetic 
radiation [I]. 

It is rather surprising that the Zitterbewegung Problem attracted only 
sporadic jSJ E] attention during years of development of the theory of quan- 
tum fields and efforts to achieve the unification of all fundamental interac- 
tions. It was demonstrated that the Zitterbewegung oscillations are 
due solely to interference between the positive- and negative- energy com- 
ponents in the wave packet. The Zitterbewegung is completely absent for 
a wave packet made up exclusively of positive energy plane wave solutions. 
It is clear from the above analysis that if one achieves the reformulation of 
the Dirac equation such that the complete orthonormal set of linearly inde- 
pendent solutions will contain only positive energy states then the Zitterbe- 
wegung oscillations will disappear. Indeed, the charge-conjugated solutions, 
associated with the positron, must be retained. 
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It has been known for a long time that the algebraic structure of Dirac 
equation is closely related to the quadratic normal division algebra of quater- 
nions. Here we suggest a quaternionic reformulation of the Dirac equations 
[S] , as well as an additional set of similar equations suitable for description of 
the free propagating quark motion. The main effort is made to obtain equa- 
tions with the intrinsic SU (2) £g> U (1) local gauge invariance. In contrast 
with the approach of S.L. Adler and others we consider the possibil- 

ity that the previously obtained quaternionic extension of the Hilbert 
space description of quantum fields represents a consistent mathematical 
framework for the electroweak unification scheme (a brief summary of rele- 
vant results is given in the Appendix). It is obvious that in order to achieve 
unification of all fundamental interactions, the algebraic extension beyond 
the quaternions is needed. We demonstrate that mathematical structure 
of the obtained equations of motion suggests that the required extension 
may proceed through wave functions which possess three and seven phases, 
whereas the scalar product remains complex. In that case the examples of 
nonextendability to octonionic quantum mechanics jU] are not valid. 

2 Equations of motion for fundamental fermions 

Let us consider the algebraic structure of the Dirac equation. The problem 
is to achieve factorization of the energy-momentum relation 



in such a way, that the correspondent Hamiltonian is the generator of 
Abelian translations in time, which is expressed by the Schrodinger equation 



It was demonstrated by P.A.M. Dirac [Q that in terms of the two- 
dimensional commutative quadratic division algebra of complex numbers 
no solution can be found. The problem requires intrinsically an extension 
of the algebraic basis of the theory. The Dirac's solution of the problem, 



77.2 2 2, 2 4 

h = p c + m c 





Hip. 





(3) 
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uses the generators of the C4 Clifford algebra: 



ajQfc + afcOj = 25ik 

oti(3 + pat = i,k = 1,2,3 (4) 

Q? = /? 2 = 1. 

However, such a drastic growth in algebra is only apparent. The true 
physical content of the obtained result is expressed more distinctly if © is 
written in the following form: 



1 dip . dip mc , 

V 

j = 1,2,3 (5) 



c dt aj dxj h 



. 1 a<p ^ . cHp mc %p 

c dt 3 dxj h 



(we choose to work in the Weyl representation |12j 
















-1 


-1 





i> \ 





P = [ i n * = 1,2,3 (6) 



by the reason, which will be explained below). 

In a precise analogy with the physical content of Maxwell's equations, 
we again have to deal with two mutually connected waves, which propagate 
together in a space with constant velocity. 

Now, the full basis of the certain algebra (C2 Clifford algebra) is symmet- 
rically used in However, it is assumed that i in © commutes with <jj, 
that is the algebra is defined over the field of complex numbers. Therefore, 
the algebraic foundation of this formulation is based on an eight-dimensional 
non-division algebra. In addition, it is customary in the applications to con- 
tinue working with complex numbers as an abstract algebra, but for the C2 
Clifford algebra, one makes use of a representation (Pauli matrices) intro- 
ducing into the theory an asymmetry, which has neither mathematical nor 
physical justification. 

Now, I will demonstrate that the algebraic foundation of (J5J) may be 
reduced to a four-dimensional real quadratic division algebra of quaternions 
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and that the structure of Dirac equations is intrinsically connected with the 
functional-analytical structures mentioned above. 
First of all, let us substitute 



e j = ~ ia j 3 = I, 2 , 3 (7) 

into ©• Then 

. 1 dtp dtp mc (p 
c dt dxj H 

J = 1,2,3 (8) 

1 d(p d(p mc 



* c dt ej dx-i h 



or, equivalently, 



e \ 1 d_ _ ( ej \ _d \ f tP 
e J c dt \ -ej J dxj ) \ tpi 



_ mc ( eo 
"IT V -e 

e \ 1 d_ _ ( ej \ d 
e J c dt ~ V - ej ) dxj 



and 



mc / e \ / tp 
IT 1 -eo M ^ 



-e A 1 /3 _ / ej \ d \ ( tP 
e ) c dt \ ej J d Xj J \ 



__ mc / e 

"I ^ e 

-eo \1 _ / ej \ d 

eo y c dt \ e-,- y <9xj 



mc / e W V 
' & V e A 



J = 1,2,3 (9) 



J = 1,2,3 (10) 
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Notice that the states have the form (Al) and all operators have the 
form (All) and (A12). 

Besides that, an additional (and only one) mass term is allowed: 





eo 







ld_ 

cdi 

eo 






-eo 
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f 




\ - — — ( 




V e 




1 cdt \ 







+ 



d 

dxj 



e 
eo 



_d_ 

dxj 



1,2,3 



eo 






-eo 



+ 



7712 C 



e 
eo 



\ V e y c a* V e i / dxj j V / 




i = 1,2,3 



(11) 

It may be verified that the energy-momentum relation is not spoiled if 
one defines 



M= Jm\ + 777^. (12) 

Here we are forced to consider masses as given phenomenological param- 
eters. If mi 7^ 0, the presence of this additional term does not increase the 
number of fundamental plane wave solutions of the equations Hll[) . There- 
fore, we will consider the equations 1)1 ljl with m\ = as a separate inde- 
pendent set and in order to maintain the direct connection with the Dirac 
equations, will neglect the 7772 term in the presence of the non-vanishing m\ 
term. 
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Indeed, only two equations (jllj) are independent: 



1 dtp . dip mc 

3 



c dt % 6,7 Oxa h 



3 = 1,2,3 (13) 
1 dcp . dcp mC ip 

c dt 1 dxj h 

The form (|13|) is very convenient for the investigation of gauge invariance 
group of the Dirac equations. The U (1) gauge invariance group from the 
right is generated by the transformations 

Ip' = TpZ, <p' = <p Z 

(14) 

z = a + bi, \z\ = 1, a, b are real numbers. 

Since for every pair of solutions (tp, (p) of the linear differential equations, 
(tpa, (pa) (a is a real number) is also a solution, it is always enough to show 
that the particular transformation 

ip = ip i 

(i.e. a = 0, b = 1) (15) 

<p' = (pi 

leaves the equations invariant. 

Invariance of the equations (|13|) with respect to this transformation is 
obvious. Let us consider what is a left gauge invariance group of the Dirac 
equations. Remember that (f 13ft are the equations for free propagating waves 
and thus admit solutions of the form 



-i(Et - px) 

w = Ul exp 

n 

-i (Et — px) 



(16) 



U2 exp ■ 



h 

Therefore, the U (1) transformations 

ip' = ziip, <p' = z x 4> 

z\ = a + bii, \zi\ = 1 

eipi + e 2 p2 + eie 2 p 3 

11 s m 



(17) 
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leave the equations (|13j) invariant and constitute the left gauge invariance 
group of the Dirac equations. 

In order to see that, it is sufficient to show again that 



ip = i\ip 

4> = h4> 



(18) 



is a solution of the equations (|13|) : 



c 

hU 2 



- (—iE) i - e\i\U\ (ipi) - e 2 hUi (ip 2 ) - e\c 2 i\\J\ (ip 3 ) = mc i\U 2 

(iE) i - e.\i\\J 2 (ipi) - e 2 i\U 2 (ip 2 ) - e\e 2 i\XJ 2 (ip 3 ) = -mc i\U\. 

(19) 



Then 



i x U x E 
c 

hU 2 E 



{eipi + e 2 p 2 + eie 2 ps) i\JJ\ i = mc i\U 2 



(eipi + e 2 p 2 + e\e 2 p 3 ) i\U 2 i = -mc i\U\. 



(20) 



By definition (see (jl7|) ). 

eipi + e 2 p 2 + eie 2 p 3 = i\ |p| . 



(21) 



Therefore, 



UiE 



U 2 E 



(eipi + e 2 p 2 + eie 2 p 3 ) XJ\i 



(eipi + e 2 p 2 + eie 2 p 3 ) U 2 i 



i\ (mcU 2 ) 



i\ {—mc U\) . 



(22) 



It is assumed in the Dirac equations ^3], that 



0, j = 1,2,3 



(23) 



and hence the obtained gauge invariance group is U (1) ®U (1). Now it be- 
comes clear why the Dirac equations allow us to incorporate an additional 
charge |14j and turn out to be suitable for the realization of the electroweak 
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unification scheme without contradiction with the Aharonov-Bohm ef- 
fect PS]. 

However, the group-theoretical content of this scheme |15j . side by side 
with the functional-analytical structures suggests that the left gauge 
invariance group should be larger (at least U (l;q) = SU (2)) and should 
not contain an Abelian invariant subgroup. The simplest way to satisfy 
these requirements is to identify the Abelian groups (|14jl and (fTTj) discussed 
above, that is to attach to the Dirac equations the following form: 



1 dif) dip mc 



c dt " dxj 3 h 



J = 1,2,3 (24) 
1 d(f> . d(f> mC ip 
c dt dxj 3 h 

and to drop the assumption (|23|) . Then the algebraic foundation of the 
theory is reduced to a four-dimensional real quadratic division algebra of 
quaternions. 

The U (1) right gauge invariance of the equations (|24|) may be main- 
tained if 



eipi + e 2 p2 + eie 2 P3 



(25) 



\P\ 

and may be demonstrated exactly in the same way as (|18|) - (|22j) . 

Consequently, we have obtained the additional meaning for i, which ap- 
pears originally in the Schrodinger equation. An algebra itself forms a vector 
space, and the basis of algebra constitutes a suitable set of orthogonal axes 
in that space, for example, a complex algebra may be considered as a two- 
dimensional plane with the orthogonal directions 1 and i. In that space 
i standing in the left-hand side of the Schrodinger equation define the di- 
rection of the time translations, which form an Abelian group. Therefore, 
the £7(1) right gauge invariance of the equations ()24|) leads us to the con- 
clusion that these equations define (the condition (|25ft ) the direction of the 
time translations at the three-dimensional quaternionic surface (the space 
of quantum mechanical phases). Then, a possible physical interpretation is 
that, compared with a classical relativistic particle, a quantum particle has 
not only its proper time but, in addition, a proper direction of time. Per- 
haps, this may serve as an explanation of why quantum equations of motion 
contain the first time derivative, whereas the classical equations of motion 
are expressed in terms of the second derivative. 
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We investigate now how our manipulations have affected the correspond- 
ing solutions. As it is well known, the general solution of Dirac equation 
may be formed as a linear combination of the four independent solutions, 
which are four spinors with four components. Two of them are obtained for 

E > for two spin states U® = ( M and V® = ( J ) , respectively. 

The other two we are forced to obtain using E < since there are no other 

possibilities. They correspond to the arbitrary choice of = f J ^ and 

Let us check what happens in our quaternionic version of the Dirac 
equation. In order to maintain connection with the original Dirac solutions, 
let us form a complete orthonormal set of it: 

M _ ( Ur (p) \ -i (Et - px) 



Ui (p)i J ™* h 

= ( u* (£ ) exp -a (Et - m 



U 2 (p)i J ^ h 
m ( U X (p) \ -i (Et - px) 



( ' ;> = ( -Ui (p)i) exp n (26) 

M _ ( U 2 (p) \ -i (Et - px) 

- V -U2 (p)i ) p n 



we have 



+ \p\ ) U x - mc U 2 = 



E 

c 

-mcUx + - \p\\ U 2 = 0. 
The existence of non-trivial solutions is ensured by 

E? 
c 2 

and 



(27) 



\p\ 2 — m 2 c 2 = 



U?> 2) = — ^E^. (28) 
1 E + c\p\ 2 K > 
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Let 



Then 



U. 



ill 



and U. 



(.2.) 



m _ mc 2 iVi ( 1 

r£) 



mc 2 Ni 



1 \ —i (Et — px) 
exp 



E + c\p\ 



i 


V o y 



cxp ■ 



-i {Et — px) 

h 



cxp ■ 



-i (Et — px) 

h 





1 

V i J 



cxp ■ 



-i (Et — px) 

h 



Now 



Let 



^(3,4) = (E + c\p\) u{3A) ^ 



U 



i3.i 



and U 



Then 



iV 2 





>( 


1 


/ 


1 


\ 




— i 











V 





/ 



exp 



-i (Et — px) 



cxp • 



-i (Et — px) 

h 



M _ (E + c\p\)N 2 ( 



mc* 



(E + c\p\)N 2 



1 / \ -i 
( \ 



1 \ —i (Et — px) 
exp — 




1 

V -ij 



exp • 



-i (Et — px) 

h 
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Our choice is made in order to compare with the standard set of the lin- 
early independent solutions for the Dirac equation. Indeed, the alternative 



*9 



mc 2 Ni 
E + c\p\ 



( 1 \ 

— i 



V o / 



exp ■ 



-i (Et — px) 

h 



and 



/ \ 



1 

\-iJ 



(32) 



exp ■ 



-i (Et — px) 



may serve us equally well and at the same time make things more trans- 
parent, since we have obtained exactly the same solutions as4» and ^f, 
which would be negative energy solutions in the Dirac equation with p'= -p 
if we make the substitution 



gipi + e 2 p2 + eie 2 P3 
\P\ 



Obviously, the obtained set is mutually orthogonal. 
Finally, using standard normalization condition, we obtain: 



(33) 



1 



2 V E + c\p\ 



1 E + c\p\ 

2 V mc 2 



1 



2 V E + c\p\ 



1 E + c\p\ 

2 V mc 2 



( 1 \ 

i 



V o / 
f°\ 



1 

V i J 

I 1 \ 

— % 

o 

V o ) 
( \ 



1 



exp 



exp 



-i (Et — px) 

h 



-i (Et — px) 



exp • 



exp 



-i (Et — px) 

h 



-i (Et — px) 



(34) 
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The obtained solutions maintain symmetry with respect to space coordi- 
nates that may be expected based on the assumption of homogeneity of the 
space-time continuum. Indeed, the correctness of the suggested equations 
may be verified only through careful comparison with the experimental data. 

Now let us consider similar equations 



1 dip . dip mc 

cp i 

J = 1,2,3 (35) 



c dt 1 dxj ^ h 



1 d(j) . d(p mc . 

c dt dxj h 

and verify that they admit an additional set of plane wave solutions, for 
example, in the following form: 

-ej (Et - px) 
ipj = Uiejexp — J - 

3 = 1,2,3 (36) 



-ej (Et — px) 
' '^^ h 



U2 exp • 



Here U\ and U2 are assumed to be real numbers. Then 



or 



/ e«Ei \ 
Uiej I — ^— - ejpiei - ejp 2 e2 -ejP 3 eie 2 l = mcU 2 i 

U 2 ~ e iV\£\ - ejp 2 e2 -ejp 3 eie 2 l = mcLq eji 



l.'\Cj [ _£z^l _ e . (eipx + C2P2 + eie2/'3 ! ) ~~ /?'<•( hi 
' e-Ei 

U 2 ( — — - ej (eipi + e 2 p2 + ei^Pz) ) = mcU\ eji. 



(37) 



(38) 



But according to (|25|) 

eipi + e 2 p2 + eie 2 P3 = i \p\ ■ (39) 
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which gives 



Ui + \p\*j i = mcU2i 

( E \ 
U2&j ( — — \p\ I i = mcUieji 



(40) 



or 



r, ( | + \p\ 



2<j \ ^ - \P\ 



and 




mcU2 
= mcUiej 

mcl/2 

= mcUiej. 



(41) 



(42) 



Thus, we finally obtain 

' E 



+ \p\ ) U\ - mcl/2 



E 



mcUx + \p\ ) U 2 = 



(43) 



which justifies the above-made assumption concerning the reality of U\ and 
U2- The existence of non-trivial solutions is ensured by 



E 2 



2„2 



— \p\ — m c 







and 



mc 



U 2 



(44) 



E + c \p\ 

Thus, we have obtained a triplet of solutions, each one associated with 



the same mass, but with 



fy, <f> k ] + j,k = 1,2,3 
V'j, ^ j / fe 



(45) 
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In addition, in the capacity of Cj, j = 1,2,3, one may choose not only 
the quaternion basis itself but other sets, for example, 

ie\i, ie 2 i, ie 2 eii (46) 

or 

eiiei, e 2 ze 2 , eie 2 zeie 2 . (47) 
Notice, however, that (|lTj) do not form a quaternion but 

eiiei +e 2 ie 2 + eie 2 ieie 2 = i. (48) 

Additional knowledge is required in order to define which set is relevant 
and how it may be associated with the correspondent physical objects. 

So far, our discussion has been restricted to the four-dimensional real 
quadratic division algebra of quaternions. However, it is clear that the 
equations ()24j) and ()35j) serve in an uniform manner also the octonionic 
extension of the complex Hilbert space . 

If the underlying algeraic foundation of the theory is extended to include 

the eight-dimensional real quadratic division algebra of octonions, then 
the corresponding additional set of solutions for the equations (|3l)|) may be 
obtained. 

They may have, e.g., the following form: 



-i (Et -px) 
tp = V\i exp — 



U2 exp 



h 

-i (Et — px) 

h 



(i is given by (|23|) 1 



, TT -ejfe {Et - px) 

Vk = Uie k exp 

n 

, TT -efc (Et - px) 
4> k = C/ 2 exp 

• % + k = 4,5,6,7 49 

, TT ■ ~3k (Et - px) K J 

ipj k = Utjk exp 

, rr ~3k (Et - px) 
4> jk = U 2 exp 

where U\, J7 2 
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are real numbers; 



3k = e k i, k = 4,5,6,7 (50) 

and, as before, form a quaternion. This quaternion turns out to form an 
algebraic foundation of the momentum space and, therefore, the algebraic 
symmetry between coordinate and momentum spaces may be broken in this 
formulation. 

The particularly symmetric case occurs, if k = 7. Then 

eipi + e 2 p2 + e 3 p 3 
\P\ 

(51) 

e^Pi + &5P2 + e 6 p 3 
3 \P\ 

Indeed, in each case the set of equations (|2*i)) and (|3l))) should be supple- 
mented by corresponding leptonic equations, for example, for (|51|) 



1 d\ . 


dx 


dx 


dx 


c dt J 


dx 


oy 


dz 


1 d£ . 


- ^e 4 - 


d£ 




~cdi J 


dx 


oy 


— e 6 

dz 



mc 



(52) 



mc 



Based on the results of M. Zorn that each automorphism of the 
octonion algebra is completely defined by the images of three "independent" 
basis units JH], it was demonstrated by M. Giinaydin and F. Giirsey ^2] 
that under given automorphism a we have three quaternionic planes in the 
space formed by octonion algebra (space of quantum mechanical phases), 
which undergo rotations by the angles 4> 1 , (f> 2 , 4> 3 , respectively, such that 

4>i + 4> 2 + </>3 = mod 2-7T (53) 

remains invariant. The planes (ej, Cj) are determined by the conditions 
Ciej = et and e k is the fixed point common to all of those planes (com- 
pare with Q49[) and ()5UjO . These results might help to extract the set of 
independent solutions and to obtain its correct classification. 
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It is worth mentioning that an alternative arrangement can be also pos- 
sible. We may consider a septet of solutions, each one associated with the 
same mass, namely 



ip = £7ie fc exp 



-e k (Et - px) 

<p = U2 exp 

n 

and, perhaps, an additional one in the form 



1,...,7 (54) 



ip = Uiejexp 

-ej (Et — px) 



ej (Et — px) 

11 j = l,-,7 (55) 



U2 exp ■ 



where 



3k = e k ie k 

jk+3 = e k+ 3 jek+3 k = 1,2,3 (56) 

37 = e 7 



or 



Jk = 

jk+3 = jek+3 j k = 1,2,3 (57) 

37 = e 7 

with the common fixed point ej: i and j are given by (|51|) . 

The number of independent solutions, which are arranged in such a way, 
is sharply reduced and serves as a slight reminder of a similar possibility 
discussed in the literature 20 . 

Now it may be clarified why we have chosen to discuss the Dirac equa- 
tions in the Weyl representation. The reason is merely technical. In order 
to perform clean octonionic calculations, we assume that the solutions of 
the equations ()35|) have the form (J36|). where U\ and U2 are real numbers. 
Then the obtained relations (|43|) justify the assumption. Notice that the 
electroweak unification scheme |15| is based on the use of this representa- 
tion of the Dirac equations; in addition, in that case the solutions behave 
naturally with respect to the Lorentz transformations. 

Let us demonstrate that the Dirac equations in the form (|24|) . as well 
as the set (|35(l . permit a consistent probabilistic interpretation (here the 
discussion is restricted to the case where the underlying algebraic basis are 
quaternions) . 
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Consider formally 



1 dtp dtp m\c m 2 c 



c * l ~ dx~^ = + 



Then 



1 dtp dtp . ra\c ± . i 7«2C 



»4 1 a ' e i l = ^ ( f )l + 

c at ax, n 



1 d<p dcp . m-\_c m2C tp 

c dt dxj 3 h h 



and 



1 dtp . dtp m\c - m 2 c- 

3 



c dt ie ° dx-i H h 



j = 1,2,3 (58) 



J = 1,2,3 (59) 



j = 1,2,3 (60) 



1 30 . dcp m\c .- m 2 c- 



c dt 3 dxj h h ^ 



Thus, we have 



1 dtp 



dtp 



c dt dxj 



mic - m 2 c - 
■—(pitp -\ —<ptp 



1 dtp 



dtp 



- ^TST + ^ dx~ = 



c ^ dt 
1 dcp- d(p - 

7 W° - dx~ e ^ 



mic - m 2 c — 
-tpicp + ——tp(p 



h 

vn\c 



tp icp 



h 

m 2 c 
IT 



j = 1,2,3 (61) 



tp<p 



1 dj> ,. d^ 
c % " dx~ 



— W Y 
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Adding all equations (|6*T]). we find 



1 f dip - , (9-0 50 -7 ,90 

- < TT-^ + '0-^- + TT-0 + 077^ 

dt* * M dV v dt 



c 



3 = 1,2,3 



(62) 

Notice that the mass terms vanish separately and, hence, the derivation 
holds separately for (|24|) and (|35l) and from now on it is understood that we 
consider the solutions of these equations separately. 

The equation Q62J1 is invariant under the gauge transformation (|15jl with 
% given by (|25[). 

Then we have 

1 9 ,7 1 5 ,7 

- 777 W + - 777 W 

C OT c <9i 

J = 1,2,3 

[ dip . - dip d(j) - dcj)\ 

+ < 77— «W + 77— - 77—^0 - 0ejZ 7— > = 0. 

(63) 

Adding equations (|6*2*)l and (|E3*|) . we obtain 

hit w + w 

1 r dip , . . . -. , , . . . dip 
+ 2^dx~^ j% + %ej > ^ + + ' 6j) J = 1,2,3 (64) 

<90 - 30 

or 

where 

p = ip^p + 00 



^ + divj = (65) 



c — — 

3k = 7| { "0 (eki + ie fc ) tp - 4> (e k i + ie k ) 0} . 



lb = 1,2, 3 (66) 
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If i is given by (|25|). then 

e k i + ie k = k= 1,2,3. (67) 

|PI 

Thus, in accordance with the Galileo, Maxwell and Schrodinger theories, 
the probability current j is proportional to the velocity operator, which is a 
constant of motion for free particles. 

Zitterbewegung [3] phenomenon is absent in this formulation. 



3 Conclusion 

In account of the experimental information that became available during 
the 

last century, it is desirable for the equations of motion for the fundamen- 
tal fermions to have the following properties: 

- the equations should possess SU (2) g) U (1) local gauge invariance 
intrinsically since the electron is not a source of pure electromagnetic 
radiation but also has the ability to participate in weak interactions; 

- the electron is the only particular member of the entire family of fun- 
damental fermions, it is desirable that all fermions are described in the 
uniform manner; 

- there exist three replication of the families of the fundamental fermions; 

- leptons do not have a color; 

- quarks do have a color; 

- each quark appear in triplet associated with the same mass; 

- a color is associated with the internal degree of freedom; 

- a color symmetry can't be broken. 

In contrast to the existing quaternionic formulations [9, 10] of the Dirac 
equation, we suggested here closely related but essentially different sets 
of equations that allow description of the free motion for electron and 
neutrinos, as well as triplets of quarks. The suggested solution possesses 
SU (2) (g) U (1) local gauge invariance intrinsically. All obtained solutions 
have the structure (Al) and (A2). These solutions are substantially different 
from the standard ones and may yield, therefore, different values for observ- 
able physical quantities. We propose to compare and verify them against the 
existing experimental information. The suggested reexamination may help 
to decide what is a relevant mathematical framework suitable to achieve a 
solution of the unification problem for the fundamental interactions. 
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During the last decades, an enormous progress in the understanding 
of quantum theory of fields took place. It became almost apparent that 
we have to deal with four essentially fundamental interactions, which have 
a similar origin, namely, presence of phases in the quantum mechanical 
description of the fundamental sources of these fields [21]. In addition, the 
investigation of the properties of the fundamental sources (leptons, quarks) 
clearly established that the quantum numbers in addition to electric charge 
(week hypercharge and color) appear in amazing correspondence with the 
complex variety of the radiated fields. 

The formulation of classical mechanics and the classical theory of fields 
have demonstrated that the presence of additional interactions requires a 
suitable generalization of the mathematical language used. Application of 
the analogy with the structure of classical physics in the framework of func- 
tional analysis naturally concentrated around attempts to extent it on all 
Hurwitz algebras, as the underlying algebraic foundation of the theory. 

The above discussion may be considered as an additional step towards 
realization of a program initiated by E. Schrodinger [22] to treat all of the 
physics as wave mechanics: 

a) The universal mathematical architecture of the physics is given in 
terms of ten functional - analytical frameworks, suitable to incorporate the 
results of physical measurements. 

Real, complex, quaternion and octonion states with real scalar product 
should be equivalent to the theory of classical fields. Unification of elec- 
tromagnetism with gravitation should occur already in the classical field 
theory. 

Complex, quaternion and octonion states with complex scalar product 
should allow realization of present unification schemes. Notice that pure 
relativistic quantum electrodynamics does not exist because there are no 
elementary sources of pure electromagnetic radiation. Neutrino is an ele- 
mentary source of pure weak radiation. 

Quaternion and octonion states with quaternion scalar product should 
describe wave mechanics of space-time continuum. 

Octonion states with octonion scalar product should allow ultimate re- 
alization of idea of elementary particles picture of natural phenomena. 

b) One expects that the quantum mechanical space-time continuum 
should be different from its classical counterpart. Perhaps, the spin is not 
a dynamical variable, but the feature of the quantum mechanical world, 
namely, the world point is described by the following expression (before 
inclusion of quantum gravity): 
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X = 



t 



-eix - e 2 y - e 3 z 
t 



e±x + e 2 y + e 3 z 



It is interesting to construct the correspondent metric space and the 
application of Least Action should lead to the equations of motion for the 
fundamental fermions. 

c) In the present discussion we consider masses as external phenomeno- 
logical parameters. However, the structure of octonion quantum mechanics 
with complex scalar product suggest natural mechanism to generate masses 
of the fundamental fermions as energy gaps obtained after splitting states 
of initially degenerated two-level physical system. 

I am grateful to L.P. Horwitz, Y. Aharonov, S. Nussinov, and I.D. Vagner 
for the stimulating discussions. 

This work has been supported in part by the Binational Science Foun- 
dation, Jerusalem. 

Appendix 

This paper is concerned with the relativistic dynamics of single particle 
states and for this reason we have used only the following results of the 
particular realization of this program for the quaternionic and octonionic 
Hilbert spaces with complex scalar products: 

1) In the quaternionic extension, quantum mechanical states are repre- 
sented by 



where f = fo + fi e i'i fo, fi are re& l functions of the space-time 

i = 1 

coordinates and ej, i = 1, 2, 3 form a basis for the real quadratic division 
algebra of quaternions; 




or 



(Al) 





3 
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In the octonionic extension quantum mechanical states are represented by: 



or 



(3) 



(4) _ 



— ( f ) 



(A2) 



V2 



i = 1 

i = 1, .. 



/o, fi are real functions of the space-time 
, 7 are a basis for the real quadratic division 



where / = / + 

coordinates and ei 
algebra of octonions. 

The ei and ej in the definition of the states (Al) and (A2) play the role 
of a label for the generator of a complex field in the space of one-body states. 
For example, any one of the quaternionic units or some linear combination 
of them 

aei + be 2 + ce 3 



v 7 ? 



(A3) 



+ b 2 + c 2 

(a, b, c are arbitrary real numbers) may be used for this purpose. Thus, a 
definition of this combination cannot be obtained kinematically and turns 
out to be a matter of the dynamics of single particle states. 

2) Consider the general form of operators, induced by the structure (Al) 
of the vector space. 

For the complex linear operators 



A, 



a n a 12 

0.21 Q-22 



(A4) 



where matrix elements are real operators over quaternions and, in turn, 
are assumed to be at least z-linear operators, we have 



^'-tAD-^-tAz ::)(/!) (A5) 

and 

«2i/ + cWiei = /'ei = (a n f + a 12 fei) a = a n fei - a 12 f. (A6) 
Therefore, 

a\2 = -«2i; an = a 2 2- (A7) 
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The restrictions (A7) on the matrix elements of the operator (A4), ob 

ned for the s 
transformations 



taincd lor the stales of the form \I>jP , are also valid if one considers the 



*f'=A z *f. (A8) 

However, for the operators transforming the state into the state 
(and vice versa), 

•r-M £ )- ( au au ) M/ )■ < A£ » 

we have 

021 = ai2; an = -022- (A10) 

Thus, we have obtained two possible types of complex linear operators, 
either 

A m = ( ou «i2 ^ (A11) 

or 

# = ( au ° 12 1 (A12) 
We remark that the matrix elements an and ai2 do not commute. 
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